Ferromagnetism of transition metals and screened exchange interactions 
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We have investigated the magnetic properties of Fe, Co, and Ni at finite temperatures on the 
basis of the first-principles dynamical coherent potential approximation (CPA) in order to clarify 
the role of the exchange interaction energy (J) screened by sp electrons and its applicability to 
finite-temperature magnetism. With use of the atomic J, we obtained the Curie temperatures (Tc) 
1930 K for Fe and 2550 K for fee Co, which are overestimated by a factor of 1.8 as compared 
with the experimental values, while we obtained Tc = 620 K for Ni being in good agreement 
with the experiment. Calculated effective Bohr magneton numbers also quantitatively agree with 
the experimental values. By comparing the results with those obtained by the screened J and by 
comparing them with the experiments, we found that the screened J, which are reduced by 30% as 
compared with the atomic ones, improve the ground-state magnetizations and densities of states at 
low temperatures in Fe and fee Co, as well as Tc in fee Co. But the screened J yield worse results 
for finite temperature properties of Fe and underestimate both the ground-state magnetization and 
Tc in case of Ni. We discuss possible origins for these inconsistencies. 



INTRODUCTION 

The band theory of itinerant magnetism has been much 
developed in the past half a century on the basis of the 
density functional theory (DFT) [lj-yl with local spin 
density approximation (LSDA) [4| or the generalized gra- 
dient approximation (GGA) [5j. The DFT explains quan- 
titatively the ground-state magnetism of 3d transition 
metals and alloys. The first-principles band calculations 
with use of the LSDA, for example, yield the magne- 
tization per atom 2.15 /xr for Fe and 0.59 /Ur for Ni, 
which are in good agreement with the experimental val- 
ues [10 2.22 /i B and 0.62 /x B , respectively. The LSDA 
stabilizes artificially the fee Fe instead of the bec Fe with 
the volume contraction by 15 %. The GGA has solved 
the problem on the stability of structure in Fe @|. 

On the other hand, understanding of the finite- 
temperature magnetism in transition metals is still far 
from the final goal of quantitative description. The 
Stoner theory based on the DFT band calculations yields 
the Curie temperatures (To), 6000 K for Fe and 3000 K 
for Ni dEil, which are 6 or 5 times as large as the ob- 
served values (1040 K for Fe (HI and 630 K for Ni [Hj]). 
Because of a large discrepancy between the theory and 
the experiment, many theories which take into account 
spin fluctuations at finite temperatures have been devel- 
oped 13]. Hubbard 14 1 and Hasegawa .15] proposed a 
single-site spin fluctuation theory (SSF) on the basis of 
the functional integral method Il6ljl9ll and the coherent 
potential approximation (CPA) [20j,|2lJ. They explained 
the Curie- Weiss susceptibility in 3d transition metals and 
obtained the Curie temperatures in Fe and Ni which qual- 
itatively agree with the experiments. 



The SSF, however, reduces to the Hartree-Fock the- 
ory at T — because it is based on a high-temperature 
approximation. Therefore, the theory does not take 
into account electron correlations at the ground state as 
discussed by Gutzwiller 22H24j. Hubbard 25|, |26|, and 
Kanamori [27]]. Kakehashi and Fulde [28| proposed a 
variational theory at finite temperatures which takes into 
account the Gutzwiller-type electron correlations at the 
ground state and reduces to the SSF at high tempera- 
tures. They found that the Curie temperatures in Fe 
and Ni are reduced by 30% ~ 50% due to correlated mo- 
tion of electrons. Hasegawa [29j also proposed a similar 
theory on the basis of the slave-boson functional inte- 
gral method. These theories are not suitable for system- 
atic improvements of the theory because of the difficulty 
in finding a suitable temperature dependence of the ef- 
fective potential being projected onto the static fields. 
In order to make a systematic improvements possible, 
we proposed the dynamical CPA [30| which fully takes 
into account the single-site spin and charge fluctuations 
self-consistently, and clarified the basic properties of the 
theory with use of a Monte-Carlo method. In the next 
paper [3l( | which we refer to I, we proposed more ana- 
lytic method to the dynamical CPA using the harmonic 
approximation (HA) |32l l33j|. 

Towards quantitative calculations, it is indispensable 
to take into account realistic band structures in solids. In 
our recent paper [34j which we refer to II, we proposed the 
first-principles dynamical CPA combining the dynamical 
CPA+HA with the LDA+U scheme (35( based on the 
tight-binding linear muffin-tin orbitals (TB-LMTO) pjij] . 
Quite recently we extended the calculations taking into 
account the dynamical corrections up to the 4th or- 
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der in Coulomb interactions, and have shown that the 
single-particle excitation spectra obtained by the first- 
principles dynamical CPA explain systematic change of 
the X-ray photoemission spectroscopy (XPS) and the 
bremsstrahlung isochromat (BIS) data from Sc to Cu at 



high-temperatures (371. 13 



In our numerical calculations, we adopted the Coulomb 
interactions U obtained by Anisimov et. al. (39| . and by 
Bandyopadhyay and Sarma [40j ] - The exchange inter- 
action energies J were taken from the values obtained 
by the Hartree-Fock atomic calculations 4l|. Recently, 



considerable efforts to obtain the first-principles U and J 
have been made at the ground state. Aryasetiawan and 
coworkers [42J obtained the Coulomb and exchange inter- 
actions for d orbitals which are screened by sp electrons, 
using the random phase approximation (RPA) and the 
LMTO method. The same type of calculations has been 
performed by using the first-principles Wannier orbitals 
by Miyake and Aryasetiawan (43|. These results indi- 
cate that irrespective of transition metal elements the 
exchange energy integral J are reduced by about 30% 
as compared with the atomic ones. Although the meth- 
ods have been applied to many other compounds [44j . 
there is no systematic investigation for the influence of 
the screened J on magnetic properties. 

We investigate in this paper the ferromagnetic proper- 
ties of Fe, Co, and Ni on the basis of the first-principle 
dynamical CPA using both the atomic exchange inter- 
action energies (atomic J) and the exchange interaction 
energies screened by 30% (screened J). We examine the 
role of the exchange energy parameters in the magnetic 
properties at finite temperatures and clarify their valid- 
ity to the magnetic properties. We will demonstrate that 
use of the screened J does not necessary improve theoret- 
ical results at finite temperatures, though we find some 
improvements at low temperature regime. 

As we have proven in our previous papers [Hi . |45[ , the 
dynamical CPA is equivalent to the many-body CPA [Ifl 
in disordered alloys, the dynamical mean-field theory 
(DMFT) [47145H in the metal-insulator transitions, and 
the projection operator CPA (EH for excitation prob- 
lem in solids. The first-principles DMFT calculations 
for Fe and Ni have been performed at the ground state 
by Miura and Fujiwara 53] within the iterative perturba- 
tion method. The finite-temperature DMFT calculations 
for Fe and Ni have been performed by usin g th e Hamil- 
tonian without transverse spin fluctuations |54j . In both 
cases, atomic exchange interaction energy J = 0.066 Ry 
was adopted for Fe and Ni. We present here the finite- 
temperature results obtained by the Hamiltonian with 
transverse spin fluctuations and by using both the atomic 
J and the screened J values. 

The outline of the paper is as follows. In the follow- 
ing section, we summarize the first-principles dynamical 
CPA based on the HA. In Sec. Ill, we present numerical 
results of calculations for the densities of states as the 



excitation spectra, magnetization vs temperature curves, 
Curie temperatures, paramagnetic susceptibilities, effec- 
tive Bohr magneton numbers, as well as the amplitudes 
of local magnetic moments in Fe, Co, and Ni. By com- 
paring these quantities calculated by the atomic J with 
those obtained by the screened J, we will clarify the role 
of the screening effects on the exchange energy J, and dis- 
cuss the availability of the screened J. Furthermore, we 
discuss the validity of the screened J by comparing calcu- 
lated magnetic properties with those in experiments. In 
Sec. IV, we present a summary of our works, and discuss 
possible reasons for the disagreement with experiments 
when the screened J are applied. 



FIRST-PRINCIPLES DYNAMICAL CPA 

In the first-principles dynamical CPA, we adopt 
the TB-LMTO Hamiltonian combined with a LDA+U 
Coulomb interactions as follows 34]. 
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We assumed here a transition metal with an atom per 
unit cell. e\ in Eq. ^ is an atomic level on site i and 
orbital L, fi is the chemical potential, i^jX' is a transfer 
integral between orbitals iL and jL' . L — (I, m) denotes 



s, p, and d orbitals. 
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ifliLc) is the creation (annihi- 



lation) operator for an electron with orbital L and spin 
a on site i, and fiiLa — a\ L<J aiLo is a charge density op- 
erator for electrons with orbital L and spin a on site i. 

The Coulomb interaction term Hi in Eq. consists 
of the on-site interactions between d electrons (I = 2). 
Uq (Ui) and J denote the intra-orbital (inter-orbital) 
Coulomb and exchange interactions, respectively. hn m 
(sum) with I — 2 is the charge (spin) density operator for 
d electrons on site i and orbital m. Note that the atomic 
level e° in Hq is not identical with the LDA atomic level 
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Here is the charge 

density at the ground state, £?lda i s a LDA functional 



to the intra-atomic Coulomb interactions }35l . 

In the dynamical CPA, we transform in the free energy 
the interacting Hamiltonian Hi into a one-body Hamil- 
tonian with dynamical potential v for time-dependent 
random charge and exchange fields, using the functional 
integral method 1^, 34| . Introducing a site-diagonal uni- 
form medium, i.e., a coherent potential £ into the poten- 
tial part, we expand the correction v — £ with respect to 
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sites in the free energy. The zeroth term in the expan- 
sion is the free energy for a uniform medium, The 
next term is an impurity contribution to the free energy. 
The dynamical CPA neglects the higher-order terms as- 
sociated with inter-site correlations. The free energy per 
atom is then given by 3l|, |34 1 



[1-1/(2Z+1)]J, J z = U /(2l+l)+J ± , and | = 
is a static field variable on a site. 

The effective potential E c h(£) in Eq. (|4]) consists of 
the static term E s t(£) and the dynamical correction term 
.Edyn(£) as follows. 



JbpA=^P]-r 1 ln J[U 
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(4) 



E eS (£) = E at (£) + E dyn (£). 



(5) 



Here j3 is the inverse temperature, J x = J y = Ji = The static term is given as 



-^Li(0)i ;i L4,(iw„)) - ^J^(,j_F L ^(iw n )F L ^(iu> n ) 



\ [-(U - 2Ui + J)J2 MO' - m - J)ni(Z) 2 + Jltl + J 2 Al 
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Here <fo Lcr (0) = w iCT (0) - T, LrT {iu) n ), and £ = £| + 
^Lcr(O) is a static potential given by Wl ct (0) = [(/To — 
2Z7i + J)n (TO (£) + (2Z7i - J)nj(£)]/2 - J^»o-/2, E iCT (iw n ) 
is the coherent potential for Matsubara frequency w n = 
(2ra + l)7r//3. The electron number nt,(£) for a given £ is 
expressed by means of an impurity Green function as 



(7) 



and = X^m^C^)- The impurity Green function 

GLa(^,i^n) has to be determined self-consistently. The 
explicit expression will be given later (see Eq. (I3"2"j) ). 

The coherent Green function FLo-(iui n ) in Eq. ([6]) is 
defined by 

F L a(iu n ) = [(i(J n -H - S(iw„)) _1 ] lL(T4LCT . (8) 

Here (H Q ) iLajL'a is the one-electron Hamiltonian ma- 
trix for the noninteracting Hamiltonian Hq, and 
(S(iw„))iLcrjL'cr = Y lLa {iui n )8 tj 8 LLI . 

The dynamical potential -Eyyn(£) in Eq. © has been 
obtained within the harmonic approximation (HA) [3l| — 
34j . It is based on an expansion of E dyn (£) with re- 
spect to the frequency mode of the dynamical potential 
VLerer' (iuiu), where w v = 2vtt/(3. The HA is the neglect 
of the mode-mode coupling terms in the expansion. We 
have then 



£dyn(£) = -/T'ln 



l + £(A,-l) 



(9) 



Here the determinant D y is a contribution from a dy- 
namical potential v^aa' with frequency w„, and the 



upper bar denotes a Gaussian average with respect to the 
dynamical charge and exchange field variables, £ m (iw ra ) 
and£ mQ (iw„) (a = x,y,z). 



The determinant D y is expressed as [34 
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Note that 1 in the above determinant denotes the 2x2 
unit matrix, a n (v, m) is a 2 x 2 matrix in the spin space, 
which are defined by 



a n {v, m) aa , 



v Lacjl {iuj v ) 



XVL<T>"cr""{-iu v )g L(T mt CT t{iu n ) , (12) 
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9Lcrcr'(iuj n ) = [(F L (iuj n ) 1 - Sv ) l ] a „i . (14) 
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Here <r a (a — x, y, z) are the Pauli spin matrices. 

A mm ' = U 5 mm , + (2U L - J)(l - £ mm '), fl£ m , = J(l - 

S mm ') (a = and B z mm , = U S mm > + J(l - 6 mm >). 

gLaa'{iw n ) is the impurity Green function in the static 
approximation, (F L (ioj n ))aa> = F La (ioj n )5oo> , and 5v is 
defined by (5va) aa i = v Laa i(0) - T, La -(iLJ n )5 (r(T /. 

The determinant D v (k ) m) defined by Eq. ([TT|) is ex- 
panded with respect to the dynamical potential as fol- 
lows. 

D v (k,m) = 1 + D^(k,m) + D^(k,m) + --- , (15) 
D^(k,m) = 22 v ai (v,m)v 11 (-v,m)--- 



CKi7i---a ri 7 Tl 



xv an (y, m)v ln (-is, m)D^j(i/, k, m) . 

(16) 

Here the subscripts on and 7, take 4 values 0, x, y, and 
z, and 



wo(i/,m) = - - i A TOW > Cm' (^)fe 



(17) 



In the calculations of the higher-order dynamical cor- 
rections [38| flj^jfs', k, to), we note that the coupling 

constants B^ m , = £? mm / = J(l - <Wi') are consid- 
erably smaller than A mm i and B* nm , because Uq and 
U\ 3> J. Thus we neglect the transverse potentials, 
v x (v,m) and v y (v,m). The approximation implies that 
a n (f, m)o-_ ff = 0. The determinant D l ,(k,m) in Eq. (|10[) 
is then written by the products of the single-spin compo- 
nents as 



D v (k,m) = D„t(k,m)D v ±(k,m). 



(22) 



Here D V(T (k, m) is defined by Eq. (fTTj) in which the 2x2 
unit matrices have been replaced by 1 (i.e., 1 x 1 unit 
matrices), and the 2x2 matrices a n (v,m) have been 
replaced by the lxl matrices a n (v,rri) aa . The latter is 
given by 



a n (v, to) 



0,z 

- E v a (v, m)v^(-v, m)h aia e na (v, to), (23) 



u Q (i/,m) = i [a = x,y,z) . 

m' 

(18) Here h 



Note that the subscript {0:7} of D^^(v,k,m) in Eq. 
(TTBl denotes a set of («i7i, • ■ ■ ,a n j n ). 

Substituting Eq. (IT51) into Eq. (jTU)) and taking the 
Gaussian average, we reach 



Edyn(£.) 
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Here each element of {l(k, to)} (k ~ 0, • • • , v — 1, m = 
1, • - - , 2i + 1) has a value of zero or positive integer. 
aj(k,m) takes one of 4 cases 0, 2, y, and z. j de- 
notes the j-th member of the (k,m) block with l(k,m) 
elements. P denotes a permutation of a set {(j,k,m)}; 
P{(j, k,m)} = {(j p , k p , to p )}. a p -i means a rearrange- 
ment of {ctj(k, to)} according to the inverse permutation 
P -1 . The coefficient C* m , in Eq. (J20J) is a Coulomb 
interaction defined by 



f<a = I -Anm' (a = 0) 

B mm > ( a = x > V' z ) 



(21) 



The frequency dependent factors _d|"^j(z/, k, to) in Eqs. 
(fl"6)) and (|20|) consist of a linear combination of 2n prod- 
ucts of the static Green functions. Their first few terms 



have been given in Appendix A of our paper II 34] 



e na (v, to) = g L cr(n - v)g La (n). 



(24) 



a(l 



), and we used a notation 



9La(n) = 9Laa(i^n) for simplicity. 

In order to reduce these summations, we make use of 
an asymptotic approximation 31, 38j |. 



e na (v, to) - q v (g L a(n - u) - 9La(n)) 



(25) 



where q v = /3/2irh i i. The approximation is justified in 
the high-frequency limit where gLa(n) is written as 



9La(n) 



1 
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iu n - €% + fj, - V L a(0) ' V 



(26) 



In the asymptotic approximation, we obtain 

n 

I)i n) Av,k,m) = Vflf 1 x Av,k,m) 

1=0 

xiM n ~° x ,(v,k,m) . (27) 

Here we wrote the subscript at the r.h.s. (right-hand- 
side) explicitly to avoid confusion. Note that the values 
of on and 7, arc limited to or z in the present approxi- 
mation. The spin-dependent quantities are given by 



D 



[^ } >,fc,m) = AW({o7})ffflW(«/,fc,m) , (28) 
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B«(i/,fc,m) = n^C^ + fc) 



f _)'-»/] r i_1 i r I - i 

EluT— vi II ~9U^ + k)\[l + —~g L ^ + k) 
i=o ' v j=-(i-») " ' q » 



(30) 



Here D^ a | (f , fc, m) = 1. ^({cry}) is the number of 
{cti7i} pairs such that aj = 7$ among the I pairs. Equa- 
tion ([30]) reduces to the result of the zeroth asymptotic 
approximation in our paper I [3ll ] when there is no orbital 
degeneracy. 

In the actual applications we make use of the exact 
form up to a certain order of expansion in D™ (see Eq. 
(|15[) ). and for higher-order terms we adopt an asymptotic 
form (|27p. In this way, we can take into account dynam- 
ical corrections systematically starting from both sides 
of the weak-interaction limit and the high-temperature 
limit. 

The coherent potential can be determined by the sta- 
tionary condition STcpa/STi = 0. This yields the dy- 
namical CPA equation as [HI 



(31) 



Here ( ) at the l.h.s. (left-hand-side) is a classical average 
taken with respect to the effective potential E c tf (£). The 
impurity Green function is given by 



EE 



SD 
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n L<y (iuj l )8Y, L<y (iuj l ) 

n v 

(32) 



The first term at the r.h.s. is the impurity Green func- 
tion in the static approximation, which is given by Eq. 
(fl4)) . The second term is the dynamical corrections, and 
KLa(iui) = 1 - F La (iujiy 2 SF La -(iuji)/ST, LlT (iuji). 

Solving the CPA equation (I5T|) self-consistently, we ob- 
tain the effective medium. The electron number on each 
orbital L is then calculated from 



(33) 



The chemical potential fi is determined from the condi- 
tion n e = Here n e denotes the conduction elec- 
tron number per atom. The magnetic moment is given 
by 



(/■/? 



P 



na 



(34) 



In particular, the I = 2 component of magnetic moment 
is expressed as 



(35) 



The above relation implies that the effective potential 
E e s(£) is a potential energy for a local magnetic moment 

In the numerical calculations, we took into account 
the dynamical corrections up to the second order (n < 
2) exactly, and the higher-order terms up to the fourth 
order within the asymptotic approximation. Summation 
with respect to v in Eqs. (TTT)1) and (|3"2"j) was taken up to 
v = 100 for n = 1 and 2, and up to v = 2 for n = 3, 4. 

When we solve the CPA equation (f31"|) . we adopted 
a decoupling approximation to the thermal average of 
impurity Green function [55j | for simplicity, i.e.. 



1 



xG La (q^/W), (Ci),iw n ).(36) 



Here we wrote the static exchange field £ as (£ z ,£^_) in 
order to make the decoupling approximation clearer. The 
approximation is correct up to the second moment (i.e., 
(£a)) an d arrows us to describe the thermal spin fluctu- 
ations in a simpler way. 

On the other hand, we adopted a diagonal approxima- 
tion [56( to the coherent Green function at the r.h.s. of 
Eq. 



F La (n) 



p™ A (e)de 



6 - Y>La{iOJn) - AeL 



(37) 



Here p^ DA (e) is the local density of states for the LDA 
band calculation, and Ae^ = (ej, — e L )5i2. The approxi- 
mation partly takes into account the effects of hybridiza- 
tion between different / blocks in the nonmagnetic state, 
but neglects the effects via spin polarization. 

The CPA equation (|31l) with use of the decoupling ap- 
proximation (1361) yields an approximate solution to the 
full CPA equation. For the calculations of the single- 
particle densities of states, we adop ted the following av- 
erage t-matrix approximation |2ll . [EtJ (ATA) after we 
solved Eq. (f3"Tj) with use of the decoupling approximation 

icsa. 



{GL<j{£,z-,i\,iUn)) - F L a(iLO n ) 
{GL<r{£z,£\>iw n ))Fl„{iu)n) 



^ A (iuj n ) = E LrT (iu n ) + 



(38) 



Here the coherent potential in the decoupling approxi- 
mation is used at the r.h.s., but the full average ( ) of the 
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FIG. 1: Spin-resolved d densities of states (DOS) of Fe at 
T/Tc = 0.52 for atomic J (dashed curve for the up-spin DOS, 
dotted curve for the down-spin DOS) and screened J (solid 
curve for the up-spin DOS, dot-dashed curve for the down- 
spin DOS). 
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FIG. 2: Total DOS of Fe for atomic J (dashed curve) and 
screened J (solid curve) at T/T c = 0.52. The XPS [H 
and BIS [6(J data at room temperature are shown by dot- 
ted curves. 



impurity Green function is taken. The ATA is a one-shot 
correction to the full CPA. 

The coherent potential T,l<t(z) on the real axis z = 
uj + ib is then calculated by using the Pade numerical 
analytic continuation method [H} . Here b is an infinites- 
imal positive number. The densities of states (DOS) as 
the single-particle excitations, Pl<t{u) are calculated from 
the relation, 

PLaiu) = --lmF Lu (z) . (39) 

7T 



NUMERICAL RESULTS 

In the numerical calculations, we adopted the lat- 
tice constants used by Andersen et. al. 36], and per- 
formed the LDA calculations with use of the Barth-Hcdin 
exchange-correlation potential to make the TB-LMTO 
Hamiltonian @. For Fe and Ni, we adopted average 
Coulomb interaction parameters U and the av erag e ex- 
change interactions J used by Anisimov et. al. (39(, and 
for Co we adopted U obtained by Bandyopadhyay et. 
al. [5(J and J obtained by the Hartree-Fock atomic cal- 
culations ; U = 0.169 Ry and J = 0.066 Ry for 
Fe, U = 0.245 Ry and J = 0.069 Ry for fee Co, and 
U = 0.221 Ry and J = 0.066 Ry for Ni. The exchange 
interaction energies mentioned above are basically atomic 
ones. As for the screened exchange interaction energy, we 
adopted 70% values [H|,|43| of these atomicj/; J = 0.046 
Ry for Fe, J = 0.048 Ry for fee Co, and J = 0.046 Ry 
for Ni. The intra-orbital Coulomb interaction Uq, inter- 
orbital Coulomb interaction U\ , and the exchange inter- 
action energy parameter J were calculated from U and 



J as U Q = U + 8 J/5, U x = U - 2 J/5, and J = J, using 
the relation Uq = Ui + 2 J. 

We present in Fig. 1 the d densities of states (DOS) 
in the ferromagnetic Fe. In the case of atomic J = 0.066 
Ry, the DOS for up-spin electrons consists of two peaks, 
the main peak due to e g electrons at U) w —0.15 Ry and 
the second peak due to t2 g electrons at cj « —0.30 Ry. 
We also find small humps at to w —0.5 Ry and 0.0 Ry. 
These humps are caused by both e g and t2 9 electrons. 
The DOS for down-spin electrons also show the structure 
consisting of the main peak at u> — 0.1 Ry mainly due 
to e g electrons and the second peak at lj ~ —0.05 Ry 
mainly due to t2 9 electrons. We find a small hump at 
v ~ —0.3 Ry due to t2 9 electrons. 

When we adopt the screened J = 0.046 Ry, the peaks 
of the up-spin DOS shift up towards the Fermi level by 
about 0.05 Ry, and those of the down-spin DOS shift 
down. It indicates the reduction of the magnetization. 
The hump at uj ~ —0.5 Ry shifts up and loses its weight, 
while the hump at ui « —0.3 Ry disappears. Accordingly 
the spectral weights of the two peaks are enhanced. The 
disappearance of the humps implies that the humps orig- 
inates in the multiplet excitations due to strong J. The 
resulting total DOS are shown in Fig. 2. The screened J 
enhances the main peak and shifts it to the Fermi level. 
The DOS obtained by the screened J seems to explain 
better the XPS data at room temperature [Hit ]. 

Solving the CPA equation (|3"T|) self-consistently at each 
temperature, we obtain the magnetization vs tempera- 
ture (ill — T) curves for Fe as shown in Fig. 3. The 
ground-state magnetizations M(0) obtained by an ex- 
trapolation of the M — T curves are 2.58 for J = 0.066 
Ry and 2.39 /i B for J = 0.046 Ry. The latter is in better 
agreement with experimental value Q 2.22 /zb, though it 
is still somewhat overestimated. The reduction of M(0) 
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FIG. 3: Magnetization vs temperature (M — T) curves, in- 
verse susceptibilities x _1 i an d amplitudes of local moments 
(m 2 ) 1 / 2 of Fe for atomic J (dashed curves) and screened J 
(solid curves). Corresponding curves in the static approxi- 
mation are drawn by dotted curves. The curve with higher 
(lower) Tc corresponds to the screened (atomic) J. Experi- 
mental data jsij ] of M are shown by +. 



with decreasing J is explained by the fact that the Hund- 
rule coupling J which parallels the spins of electrons in 
each atom competes with the disordering due to electron 
hopping and thus M(0) is not saturated. We also find 
the reduction of the amplitude of local moment (m 2 ) 1 / 2 
by 0.066 /xb when the screened J is used, as seen in Fig. 
3. 

With increasing temperature, the magnetization for 
screened J shows less temperature dependence at low 
temperatures and rapidly decreases near Tq. Resulting 
Tc is higher than that obtained by the atomic J. Basi- 
cally the Curie temperature Tc is determined by the ratio 
of the magnetic energy to the magnetic entropy. The en- 
hancement of Tc with reducing J in Fe is attributed to 
the reduction of the magnetic entropy with the collapse 
of the Hund rule coupling. Further reduction of J should 
decrease Tc because of the reduction of magnetic energy. 

Calculated susceptibilities of Fe follow the Curie- Weiss 
law as shown in Fig. 3. We find the effective Bohr mag- 
neton numbers m e g = 3.0 /ie for J = 0.066 Ry and 
TO c ff = 2.8 hb for J = 0.046 Ry, which are compared 
with the experimental value 62] 3.2 [1b- The screened J 
underestimates the effective Bohr magneton number. We 
have also calculated the effective Bohr magneton number 
for fee Fe at high temperatures (~ 2000 K) , and obtained 
TO c ff (fee) = 4.2 fiB for J = 0.066 Ry and m e ff(fcc) = 2.6 
Hb for J = 0.046 Ry after having included the volume 
effects on the susceptibilities. These values should be 
compared with the experimental value 63[ 7.0 /ib- The 
results again show that the screened J value underes- 
timates TO c ff. Present results indicate that the high- 
temperature properties of Fe are explained better by the 
atomic J, while the low-temperature ones are explained 
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FIG. 4: Total DOS of fee Co for atomic J (dashed curve) and 
screened J (solid curve) at T/T c = 0.39. The XPS [Hg and 
BIS [6^] data at room temperature for hep Co are shown by 
dotted curves. 



by the screened J. 

In the case of fee Co, the d DOS for up-spin electrons 
shift up by 0.04 Ry and those for down-spin electrons 
shift down by 0.01 Ry when the atomic J is replaced 
by the screened J. The total DOS shifts towards the 
Fermi level and the peak at u> = —0.1 Ry is considerably 
enhanced as shown in Fig. 4. 

It should be noted that the XPS [H and BIS 
data are obtained for the hep Co at room temperature, 
thus they are not able to be compared directly with the 
present results for fee Co. Nevertheless, we find the peak 
position of the BIS data is in accordance with the hump 
position in the DOS. The band width of the DOS seems 
to be larger than that expected from the XPS data. 

We present the magnetization vs temperature curves 
of fee Co in Fig. 5. The ground-state magnetizations 
obtained by an extrapolation are 1.82 /iB for the atomic 
J = 0.069 Ry and 1.72 fi B for the screened J = 0.048 
Ry, respectively; the latter obtained by the screened J is 
in better agreement with the experimental one Q (1.74 
/is)- The amplitudes of local moments decrease by about 
0.007 fiB with the reduction of the Hund-rule coupling by 
0.021 Ry. 

The magnetization for the screened J decreases in the 
same way as in the atomic J with increasing tempera- 
ture (see Fig. 5). The calculated T c is 2075 K for the 
screened J = 0.048 Ry, and is smaller than Tc = 2550 
K for atomic J — 0.069 Ry. The former (latter) is 1.5 
(1.8) times as large as the experimental value of fee Co 
(Tc = 1388 K) [65]. We also calculated the susceptibil- 
ities. The inverse susceptibilities follow the Curie- Weiss 
law but show an upward convexity to some extent. Cal- 
culated effective Bohr magneton numbers are 3.0 /is in 
both cases being in good agreement with the experimen- 
tal value 62j 3.15 /is- They hardly change due to the 
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FIG. 5: Magnetization vs temperature (M—T) curves, inverse 
susceptibilities an d amplitudes of local moments (m 2 ) 1 ^ 2 
of fee Co for atomic J (dashed curves) and screened J (solid 
curves). Corresponding curves in the static approximation 
are drawn by dotted curves. The curve with higher (lower) 
7~c corresponds to the atomic (screened) J. Experimental 
data [13] of M for fee Co are shown by open squares and 
dashed line. 
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FIG. 7: Magnetization vs temperature (M — T) curves, in- 
verse susceptibilities x _1 i an d amplitudes of local moments 
(m 2 ) 1 ' 2 of Ni for atomic J (dashed curves) and screened J 
(solid curves). Corresponding curves in the static approxi- 
mation are drawn by dotted curves. The curve with higher 
(lower) Tc corresponds to the atomic (screened) J. Experi- 
mental data [6y] of M are shown by +. 
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FIG. 6: Total DOS of Ni for atomic J (dashed curve) 
screened J (solid curve) at T/T c = 3.2. The XPS [H| 
BIS [6(| data are shown by dotted curves. 
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reduction of J. These results indicate that effects of the 
Hund-rule coupling on the local magnetic moment are 
less significant as compared with the case of Fe because 
the number of d holes is rather small (about 2.36). 

The DOS for Ni are presented in Fig. 6. They are 
shown for paramagnetic state because calculated Tq of 
Ni is lower than those in Fe and Co by a factor of 3 or 5, 
and the calculations of DOS at low temperatures below 
Tc are not easy in the present method with use of the 
numerical analytic continuation. The paramagnetic DOS 
are not sensitive to the choice of J since the Hund-rule 
coupling is not effective for the system with small number 
of d holes (w 1.32). Calculated ri-band widths seem to be 



somewhat larger than those expected from the XPS data. 
The same tendency was found also in fee Co (see Fig. 
4). The results suggest that the higher-order dynamical 
corrections are desired for the quantitative description of 
the DOS in these systems. 

The magnetization vs temperature curves for Ni are 
presented in Fig. 7. The ground-state magnetization ob- 
tained by an extrapolation is 0.64 /xb for atomic J = 
0.066 Ry and 0.46 Mb for screened J = 0.046 Ry. The 
experimental value [7j 0.62 /xb is in-between. The ampli- 
tude of local moment for J = 0.066 Ry is 1.93 /Ub, for ex- 
ample, at T = Tq. The amplitude hardly changes even if 
the screened J is adopted. This is because the Hund-rule 
coupling is negligible for the system with small number 
of d holes. 

The Curie temperatures in Ni are strongly reduced by 
the dynamical effects as shown in Fig. 7. We obtained 
T c = 620 K for the atomic J (=0.066 Ry) and 355 K for 
the screened J (= 0.046 Ry). The former is in good agree- 
ment with the observed Tq = 630 K [T^. The screened 
J underestimates the Curie temperature. The results 
indicate that the value of J is significant for the spin po- 
larization and the inter-site ferromagnetic coupling in Ni. 
In order to examine the ambiguity of the Coulomb inter- 
action energy parameters, we performed the same calcu- 
lations with use of different sets of (U,J): (0.25,0.073) 
Ry and (0.25, 0.051) Ry, where U = 0.25 Ry was ob- 
tained by Bandyopadhyay et. al. [40j and J = 0.073 Ry 
was taken from the Hartree-Fock atomic calculations [41( . 
We obtained Tc = 680 K for the former (, i.e., atomic J) 
and 325 K for the latter (, i.e., screened J); the results 
are not so sensitive to the choice of U. 

Calculated inverse susceptibilities are upwards convex, 
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being in agreement with the experimental data [62j for 
Ni. Effective Bohr magneton numbers calculated at T ~ 
2000 K arc 1.6 fiB f° r both the screened and the atomic 
J. The results quantitatively agree with the experimental 
value 62] 1.6 fiB- 



SUMMARY 

We have calculated the magnetic properties of Fe, Co, 
and Ni at finite temperatures on the basis of the first- 
principles dynamical CPA+HA, and investigated the ef- 
fects of screened exchange interaction energy J on the 
magnetic properties in these ferromagnets. 

When we apply the atomic J, we obtained from the 
magnetization vs temperature curves the Curie temper- 
atures 1930K for Fe, 2550K for fee Co, and 620K for 
Ni. Although calculated 2c in Ni is close to the exper- 
imental value 12] 630K, those in Fe and Co are over- 



estimated by a factor of 1.8 as compared with the ex- 
perimental values, 1040K (Fe) [n| and 1388K (Co) \&^, 
respectively. In the present calculations using atomic J, 
the ground state magnetizations obtained by an extrap- 
olation are 2.58 /xb (Fe), 1.82 /j,b (fee Co), and 0.64 /j,b 
(Ni). Calculated magnetizations for Fe and fee Co are 
considerably larger than the experimental values, 2.22 
fi B (Fe) [6] and 1.74 jU B (Co) Q, while that for Ni agrees 
well with the experiment (0.62 /is)- The results indicate 
that the atomic J overestimate the ferromagnetism of Fe 
and fee Co at low temperatures. We have calculated the 
paramagnetic susceptibilities. Calculated effective Bohr 
magneton numbers, 3.0 /xb (Fe), 3.0 /j,b (Co), and 1.6 /j,b 
(Ni), quantitatively explain the experimental data. 

The reduction of the exchange interaction energy by 
30% due to screening by sp electrons weakens the Hund- 
rule coupling which builds up the local magnetic moment 
on each atom. We found that it reduces the ground- 
state magnetization M(0) and the amplitudes of local 
moments in Fe and Co, especially the ground-state mag- 
netizations 2.39 (is (Fe) and 1.72 /ib (fee Co) are in good 
agreement with the experiments. The screened J reduces 
the exchange splitting between the up and down DOS. 
The reduction of the splitting shifts the main peak in 
the spin-resolved DOS to the Fermi level and enhances 
the peak in case of Fe and Co. The screened J re- 
duces the humps in the DOS below Tq which originate 
in the multiplet-type of excitations associated with J. 
The screened J also weakens the inter-site magnetic cou- 
plings. It reduces the magnetization in Ni, where the 
on-site Hund-rule coupling is not significant because of a 
small number of d holes. 

With reducing J, the magnetic entropy and the mag- 
netic energy decrease in general. We found that the 30% 
screened J enhances the Curie temperature Tc by 140 K 
in Fe because of the reduction of the magnetic entropy, 
while it reduces Tc in fee Co by 480 K because of the re- 



duction of the magnetic energy. The Curie temperatures 
of Fe and Co are overestimated by a factor of 1.5 ~ 2.0 
irrespective of the choice of J in the present calculations. 
It should be attiributed to the magnetic short range order 
which is not taken into account in the present theory. In 
the case of the susceptibilities, the screened J excessively 
reduces the effective Bohr magneton number in Fe where 
the Hund-rule coupling builds up the on-site magnetic 
moment. 

By comparing both results of atomic J and screened J 
with the experimental data, we found that the screened J 
explains better the ground-state magnetization M(0) and 
the DOS at low temperatures in case of Fe and Co. But 
it yields the underestimate of M(0) in Ni. The atomic 
J seems to be better for the explanation of M(0) in Ni. 
Concerning the finite-temperature properties, we found 
that the screened J improves the Curie temperature in 
Co, while the atomic J yields better Tc in case of Fe 
and Ni. The effective Bohr magneton numbers are quan- 
titatively explained by the atomic J. The screened J 
underestimates m e ff for both the bec and fee Fe. 

The above-mentioned results indicate that only for fee 
Co the screened J seems to improve both the low- and 
the high-temperature properties. In the case of Fe, the 
high-temperature properties seem to be explained bet- 
ter by the atomic J, while the low-temperature ones 
by the screened J. In Ni, both the low- and the 
high-temperature properties are explained better by the 
atomic J. 

Although we have to improve further the theory of dy- 
namical CPA+HA at low-temperatures in order to obtain 
more solid conclusion, the present results obtained in the 
high-temperature region should be reasonable because 
the theory is an approach from the high-temperature 
limit and the magnetic short range order is not significant 
at high temperatures. The results for the effective Bohr 
magneton numbers indicate that one has to take into ac- 
count the temperature dependence of J for the system 
like Fe where the Hund-rule coupling builds a large local 
magnetic moment but competes with the kinetic energy 
of d electrons. It is plausible that spin fluctuations at fi- 
nite temperatures break the screening of J. The present 
results also suggest that the 30% screening of J is an over- 
estimate for Ni. The range of applications of screened J 
and its validity have to be examined further in the future 
investigations with use of more advanced theories. 
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